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ABSTRACT 

In this paper we revise and complete the photometric survey of the instabihty strip 
of the southern open cluster NGC 2516 published by Antonello and Mantegazza (1986). 
No variable stars with amplitudes larger than 0'".02 were found. However by means of 
an accurate analysis based on a new statistical method two groups of small amplitude 
variables have been disentangled: one with periods < 0*^.25 (probably 5 Scuti stars) and 
one with periods > 0*^.025. The position in the HR diagram and the apparent time-scale 
may suggest that the stars of the second group belong to a recently discovered new class 
of variables, named 7 Dor variables. They certainly deserve further study. We also 
present a comparison between the results of the photometric survey and the available 
pointed ROSAT observations of this cluster. 

Subject headings: Stars: variables: delta scuti - Stars: variables: other 
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1. Introduction 

The observation of pulsating variables in homoge- 
neous samples such as stellar chistcrs is a good way 
to collect information on the effect of age, chemical 
composition and rotation on pulsation. 

The studies of open clusters in the northern sky 
made by several authors (see Slovak, 1978) have 
shown that cluster variable stars of 5 Set type, i.e. 
the typical low-instability strip pulsators, have an in- 
cidence of around 30% apparently independent of age. 

The southern open cluster NGC 2516, located at 
right ascension a(2000) = 7'* 58™. 3 and declination 
(5(2000) = —60° 52' is a young cluster believed to 
have a common origin with a Per, the Pleiades and 
IC 2606 (Eggen 1983), but imhkc these clusters it has 
an unusually high abundance of Ap stars. An exten- 
sive study by Mermillod (1981) on young open clus- 
ters identifies it as the prototype of a group, named 
the NGC 2516 cluster group, that also includes NGC 
2168, NGC 2301, NGC 3114, NGC 5460, NGC 6025, 
NGC 7243. They have almost the same age (log age 
=8.3) and are all characterized by a peculiar gap near 
the turn off point and by the presence of numerous Ap 
stars in the extreme blue region of the main sequence. 

A first survey for the detection of variable stars 
in the lower part of the instability strip was carried 
out by Antonello and Mantegazza (1986, hereinafter 
paper I) and provided evidence, besides a normal in- 
cidence of shorter period variables (P < .25 d), of a 
number of longer period variables {P > .25 d), lying 
on the cool border of the instability strip. 

In the last years the presence of variable stars 
showing small amplitudes with longer time scales than 
the typical S Set periodicities has been often reported. 
Mantegazza et al (1993) first gathered a sample of 8 
field stars and 5 cluster stars with very similar loca- 
tion in the HR diagram, on or just beyond the cool 
border of the instability strip, apparently showing a 
common behavior. They have been named 7 Dor stars 
and an updated list has been recently published by 
Krisciunas and Handler (1995). 

In the present paper we present previously unpub- 
lished observations of NGC 2516 that complete the 
survey for variability inside the instability strip of the 
cluster. The survey is complete in the sense that all 
the stars in the instability strip have been searched 
for variability. 

The absence of a physical lower limit to the pos- 



sible variability of stars implies the need of a statis- 
tical method of data analysis in order to identify the 
variable stars according to a probability criterion. A 
method was developed and described in paper I. 

In this paper we present and use a different statis- 
tical approach with respect to paper I, less sensitive 
to the photometric quality of the night. A comparison 
between the two approaches is also reported. 

2. Observations 

A new set of observations were collected at the ESO 
1 m reflector from 1987 January 4 to .January 9 in the 
Johnson B filter with the same equipment of paper 
I, i.e. a single channel photometer with a EMI 6256 
photomultiplier. In the first 4 nights 6 stars per night 
(labeled si — sq) were observed following a cycle si — 
S2 — S3 — S4 — S5 — si — S2 — S3 — S4 — S6- The cycle was 
designed in order to observe in each night and for two 
consecutive nights two of the stars already discovered 
as longer period variables (see paper I), namely stars 
S5 and Sfl. Five stars were observed in each of the 
last two nights according to the cycle si — S2 — S3 — 
S4 — S5 and in this case no previously known longer 
period variables were included. Each measurement 
consisted of 25 integrations of 1 second: the number 
of integrations was increased than in paper I in order 
to reduce scintillation noise. Each observed star was 
monitored for 7.5 hours and an average number of 85 
and 75 points per star have been collected in the last 
two nights and in the first 4 nights respectively. 

In order to avoid spurious effects due to variation 
of air transparency the analysis was based on differ- 
ential time series. Since the difference in airmass be- 
tween the stars observed in each night is very small 
(< 0.005), these spurious effects due to variable air 
transparencies are significantly reduced. 

3. Statistical Analysis 

In this run, as for the observations reported in Pa- 
per I, no comparison stars were assigned a priori: we 
then do not know which of the stars in each observ- 
ing cycle are stable or variable. On the other hand 
the expected light variations have such small ampli- 
tudes (< 0'".02) that this choice cannot be based on 
a visual inspection of the differential light curves, so 
a technique had to be developed that extract this in- 
formation from the differential time series between all 
the possible pairs of stars in the cycle. 
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In paper la statistical approach based on least 
squares analysis of differential time series was adopted. 
For the n stars in each observing night the variances 
sfj of the n(n — 1) possible combinations were com- 
puted. Under the assumption that each of these vari- 
ances is the sum of the contributions of the two com- 
ponents 
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the values of such components were estimated by solv- 
ing the system of n(n — 1) equations with n unknowns 
by means of the least squares method. 

Furthermore the quantities sf were expressed as 
= af + bf + (?, where of is the intrinsic variance of 
the star brightness, bf the variance due to the statisti- 
cal fluctuations of the counts (Poisson's variance) and 
a term which depends on the photometric quality 
of the night. The hf terms could be estimated through 
the mean number of counts collected for each star and 
subtracted, but the (? term remained undetermined 
and statistics could be performed only on the quanti- 
ties 
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Variable stars were discriminated on the basis of 
an F-test on the quantities Fi = rff /(i^„j„ . 

However the least squares based procedure may 
not be the best choice. The n(n — 1) equations, for 
instance, are made up of pairs of equations, 
and {j,i) that differ only because in the former the 
times of the i-th star have been interpolated to the 
times of the j-th star while in the latter vice versa. 
The number of real independent equations is there- 
fore n(n — l)/2, resulting in poor oversampling with 
respect to the n unknowns. Furthermore if the con- 
tribution to the variance of the sky and background 

is appreciable, e.g. of the same order of the intrin- 
sic variances, a least squares procedure can become 
unstable. 

We therefore decided to make use of a different ap- 
proach that does not involve the least squares method. 
Each differential data point in the series is defined by 



^mij{tk) = -2.5 log 
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where Ii{tk) and Ij{tk) are the counts obtained at the 
sampling instant tk, and are therefore affected by a 
statistical error ^/I^{tk) and ^JIj{tk) respectively, due 



to Poisson's statistics. Error propagation states that 
contribution of the statistical error to the variance on 
the single differential measurement is the following 



= 1.179 (- 
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that can be directly subtracted from measured vari- 
ances 



'Ni 



Ni. 



, fe=i 



(5) 

where Nij stands for the number of measurements in 
the given time-series. 

The Poisson corrected variance may be expressed 

as 
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i.e the sum of the components due to each star fluctu- 
ation plus an additional term due to sky and instru- 
mentation. Indeed the term due to sky and instru- 
mentation affects directly the differential measure- 
ment and in a first approximation does not depend 
on the stars observed. Being corrected for statistical 
noise contribution and sky background the quantities 
sf may be considered equivalent to the quantities a| 
better than those named s^ in paper I: we will there- 
fore refer hereinafter to sf as the intrinsic variances. 

The information we need to disentangle variables 
from constant stars, i.e. sf, is then distributed into 
n — 1 quantities in the form (6) and is contaminated 
by a systematic background error sf . 

To extract this information we compute the sum- 
mation over each possible differential combination of 
the i-th star 



S? = i: sh = {n- 2)sf + s| + (n - l)sf 

\j=l 

(7) 

which is a quantity related to the intrinsic variance 
sf . The term between parenthesis is constant for all 
the stars observed in the night. In order to estimate 
such a constant term we summed the quantities (7) 
over all the stars observed in the night 
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S^ = Y,Sl = 2{n-l)Y,s]+n{n-l)sl. (8) 

i=l = 1 

Through eqn (7) and (8) we then obtain 

C2 c2 1 
c2 = 2 i„2 /q^i 

' n-2 2(n-l)(n-2) 2 ^ ' 

where the only unknown on the right hand side is the 
background term s^. 

The observational procedure does not allow a pre- 
cise estimate of the background term s^. However the 
minimum value among the differential variances ob- 
tained in the considered night put an upper limit to 
this quantity, i.e. 

< min{s1j) i = l,n j = l,n. (10) 

Indeed from equation (6) one can see that equals 
the minimum value if sf = s'j = 0. The minimum 
differential variances obtained range between 1 and 
1.5 10-5 mag^ 

By introducing this upper limit in (8) instead of 
the real unknown value of background uncertainty, 
we tend systematically to underestimate the intrinsic 
variances, the only risk being to assume as constant 
a variable star. For the purposes of this paper such a 
conservative character is welcome. 

We then obtained a scale of the background-corrected 
estimates of the intrinsic variance for each of the stars 
observed in a given night. As in Paper I we could de- 
lineate a limit between constant and variables stars 
through a statistical F-test by comparing the intrin- 
sic variances with the minimum value for the night. 
However it should be noticed that the application of 
an F-test here as well as in the context of the anal- 
ysis of Paper I is not rigorous: indeed the variances 
involved are not directly obtained from the original 
data, but derived through some manipulations. 

For this reason we preferred to assume a fiducial 
limit of two times the minimum intrinsic variance for 
the night, i.e. greater than the typical F-test value of 
Paper I. 

We report in tab. 1 the results of the analysis, while 

a reanalysis by means of the new procedure of data 
published in Paper I is reported in tab. 2. 



A precise estimate of the time-scale of variability 
for the recognized variables can not be obtained with 
such a short baseline. Indeed the average baseline 
of 0.3 d only allows a poor spectral resolution of 3 
d"^. However for the purposes of this paper we are 
only interested in distinguishing shorter time-scale 
variables (P < .25 d) and longer time-scale variables 
{P > .25 d): this information can in fact be extracted 
even from poor spectral analysis. 

We performed a least-squares frequency analysis 
by means of Vanicek's (1971) technique and labeled 
L the stars that provided evidence of reduction factor 
concentrated at low frequencies (/ < 4 d"^) and S all 
the other variables (see columns 3 and 6 of tab. 1 and 
2). 

4. Results 

The new set of observations confirmed the vari- 
ability of 7 stars already recognized as variables in 
the first run, 3 of them labeled L (C52, C69 and C93) 
and 4 of them labeled S (C84, C54, C55, C51). A 
further group of 3 stars (C106, C92 and C64) were 
pointed out to be variables of type L. 

As an example of the quality if the light curves 
obtained during the survey differential curves of the 
shorter period variable C55 and the longer period 
variable C69 are reported in fig. 1 together with the 
constant star C97. Each curve is computed with re- 
spect to C98, the star showing the minimum intrinsic 
variance according to the statistics used. 

Only a few discrepancies have been found between 
the results obtained with the old and new approach 
in the analysis of the first run, and between the re- 
sults of the two runs. The new analysis of the night 
HJD 2446115 provided the following different results 
with respect to paper I: the star C21, indicated as an 
S variable, was judged constant by the new statistics, 
while the star C62, previously considered constant, 
turned out to be an L variable. Such different re- 
sults are probably due to the non perfect atmospheric 
conditions on that particular night to which the new 
statistics is far less sensitive. 

The new observations of the stars C96 and C59 do 
not confirm the variability observed in the first run, 
while C38, constant in the first run, was classified as 
S variable in the second run. We notice that all the 
three stars provided evidence of the highest variance 
among those observed on the night in which they were 
classified as variable, so they were very likely varying 
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C54 
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V(S) 


C55 


1.163 


V(S) 


C69 


1.822 


V(L) 


C106 


2.263 


V(L) 


C38 


2.284 


V(S) 


C69 


6.469 


V(L) 


C108 


.290 


c 


C66 


.495 


c 


C77 


.336 


c 


C60 


.535 


c 


C59 


.396 


c 


C17 


.676 


c 


ClOl 


.452 


c 


C51 


1.827 


V(S) 


C92 


1.557 


V(L) 


C64 


5.983 


V(L) 



Table 1: Results of the statistical analysis with the 
new approach of the last season data. The variances 
are expressed in xlO~^ mag^. V stays for variable 
star while C for constant star according to the classi- 
fication reported in the text. L and S means long and 
short period variables respectively. 
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C38 


.320 
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.884 


V(S) 
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.265 
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.284 
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C77 


.202 


C 


C99 


.292 


c 


C45 


.209 


C 


C32 


.417 


c 


C50 


.227 
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C31 


.536 


c 


C53 


.384 


c 


C69 


1.699 


V(L) 


C52 


.704 


V(L) 



Table 2: Results of the statistical analysis with the 
new approach of the data presented in paper I. The 
variances are expressed in xlO~^ niag^. V stays for 
variable star while C for constant star according to the 
classification reported in the text. L and S means long 
and short period variables respectively. Objects of the 
first run that yielded different results with respect to 
the statistical approach of paper I are marked with a 
star 
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Fig. 1. — Differential light curves of the L-type vari- 
able C69 and the S-type variable C55 observed in the 
night of HJD 2446803. The star C98, the one with 
the minimum intrinsic variance in the night, was used 
ascomparison star. The differential light curve of the 
constant star C97 is also riported for comparison. 
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when observed for the first time. 

In the case of C96 this discrepancy could possi- 
bly be due to the fact that a longer period variable 
observed on a short baseline near a maximum or a 
minimum could easily be misinterpreted as constant. 
In the case of the short period variables C59 and 
C38, one possibility is that we are dealing with multi- 
periodic pulsators, as most 6 Scuti stars are, and the 
apparent constancy is due to a disruptive beat be- 
tween different modes. The possible variability of 
C38 was already discussed by Snowden (1975) and 
Maitzen and Hensberge (1981). 

The four newly discovered L type variables (C62, 
C106, C92 and C64), added to the group already 
pointed out in paper I (C93, C52, C96 and C69), 
make up a sample of 8 early F-type, longer time-scale 
variables. 

5. Discussion 

A sample of 44 A- and early F-typc stars has 
been observed in both runs. This sample is com- 
plete in the sense that all dwarfs and subgiants with 
0.00 < {b — y)o < 0.25 present in the cluster have been 
observed. Among these stars, six were pointed out 
to be S variables (13% of the sample, excluding the 
doubtful cases C59 and C38) and 8 to be L variables 
(18% of the sample). While the incidence of shorter 
period variables in this cluster, as already stated in 
paper I, can be considered normal, the incidence of 
longer period variables is certainly surprising. 

All the stars observed are reported in the HR di- 
agram shown in fig. 2 . The stars are plotted with 
their absolute visual magnitude Vo versus their dered- 
dened (6 — y)o color index. The Stromgren photome- 
try used as well as the evaluation of the color excess 
E{b—y) and visual magnitude My are those published 
by Snowden (1975): main sequence and borders of the 
instability strip are taken from Rodriguez et al (1994). 

With the exception of the peculiar case C 38 all dis- 
entangled shorter period variables fall within the lim- 
its of the 6 Scuti instability strip and they very likely 
belong to this class of variables. We notice however 
that no S Scuti candidate cooler than FO was found. 

Very little can be said on longer period variables 
but it is worth making some speculation on their pos- 
sible nature. The observed variability is very likely 
due to intrinsic changes in luminosity of the sources. 
Indeed a spurious nature due to an observational ef- 
fect, e.g. air transparency variations, would unlikely 
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Fig. 2. — Position in the HR diagram for all the stars 

observed in NGC 2516. Squares, filled squares and 
crosses represent L-type, S-type variables and non- 
variable stars respectively: those with the symbol in- 
cluded in an open circle have been found to be X-ray 
sources. Field 7 Dor objects are represented by tri- 
angles 
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pass the filter of the statistics used: each star is 
checked against five other selected nearby stars in the 
cluster. In addition if they were due to an observa- 
tional effect they woiild have been spread out in the 
HR diagram. They are however located in a narrow 
color range. 

They could be binaries: there is no way to rule 
out this hypothesis through our survey data. We no- 
tice however that there is again no reason why in a 
bias-free sample of stars the binaries should be con- 
centrated on the cold part of the instability strip as 
they appear. 

Could they be intrinsic (pulsating, spotted) vari- 
ables ? Their position in HR diagram and their ap- 
parent time-scale suggests they could be 7 Dor star. 

The 7 Dor stars are variables in time scales be- 
tween — 5 and with typical amplitudes of 
O^.OS — 0™.! lying in a very narrow region on or just 
beyond the cool border of the instability strip. Even 
if some controversy is still present, they have been 
recognized as the first evidence of non radial g-mode 
pulsation among F-type dwarfs. The position in HR 
diagram of the most thoroughly studied 7 Dor stars 
is reported (triangles) in fig. 2 for comparison. 

Not much of an interpretation can be worked out 
from survey data of the kind presented in this paper 
and the possible 7 Dor nature of the longer period 
variables should be regarded as guessed. 

6. X-ray data 

The Interaction between pulsation and chromo- 
spheric or coronal heating mechanisms has been re- 
cently indicated as a possible important tool to col- 
lect information both on pulsation mechanisms and 
heating. 

Known single X-ray sources among early F-type 
stars are generally related with emission from a corona 
heated by magnetic fields. Two important cases are 
reported in the literature: 47 Cas (Giidel et al., 1995) 
and 71 Tau (Stern et al 1994). The former is reported 
to be a powerful X-ray source showing pseudo-regular 
variability on a time scale of the order of 1 d. Its 
Stromgren photometry (Olsen 1983) places it in the 
7 Dor region of the HR diagram, and the time scale 
falls in the typical range of the 7 Dor optical vari- 
ability. However, recent photometric observations ex- 
clude light curve variability of 47 Cas larger than 4 
mmag (Poretti, private communication). 71 Tau is 
the strongest X-ray emitter in the Hyades; it is a 5 



Scuti pulsator with spectral type A9 and large rota- 
tional velocity. It is a binary star, but it is difficult to 
explain the intense X-ray emission only with a late- 
type companion. 

The simultaneous availability of our complete sur- 
vey and ROSAT data on NGC 2516 prompted us 
to look for relations between variables and X-ray 
sources. 

The open cluster NGC 2516 was observed with 
the Position Sensitive Proportional counter (PSPC) 
at the focus of the X ray telescope on board ROSAT. 
This satellite contains an X-ray telescope with a 2° 
field of view. The PSPC is a gas-filled proportional 
counter sensitive in the energy band 0.1-2.4 keV, with 
a spatial resolution of ~ 25" in the center of the fo- 
cal plane. A detailed description of the satellite and 
detectors can be found in Triimper (1983) and Pfef- 
fermann et al. (1986). 

The observation reported here was obtained on 
1992 April 6 with a total effective exposure time of 
9284 s. Due to Earth occultation, radiation belt pas- 
sages and observations of other targets, the data were 
spread over two days. The data were analyzed with 
the XANADU package: first the event files were taken 
from the ROSAT on-line archive at the Max Planck 
Institute (Munich) and images were extracted. The 
exposure maps were linearly interpolated and rescaled 
so as to be overlaid on the images. A boresight cor- 
rection of ^ 10" was applied. Exposure corrected 
count rates or 3 a upper limits were finally derived 
(see tab. 3). 

The analysis showed that none of the suspected 
variables in the cluster has an X-ray emission (0.1- 
2.4 keV) above ~ 5 x 10^^ erg s~^ (assuming a count 
rate to flux conversion factor of ICR = 10~^^ erg s"-^ 
cm~^ and a cluster distance of 440 pc). This upper 
limit has to be compared with ~ 10^^ erg for 47 
Cas (Giidel et al. 1995) and - 3 x lO^o erg s"! for 71 
Tau (Stern et al 1992). It is interesting to note that, 
despite the non-detection of variable stars, some stars 
inside the variability strip have been revealed with an 
X-ray luminosity of a few 10^^ erg s~^. 

7. Conclusions 

In this paper we have reported on a survey for vari- 
ability of a selected sample of A- and early F-type 
stars in the southern open cluster NGC 2516. The 
sample is complete in the sense that all dwarfs and 
subgiants with 0.00 < {b — y)o < 0.25 belonging to 
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Star 


Count rate or 3 cr upper limit 


S/N 


Star 


Count rate or 3 ct upper limit 


S/N 


C4 


(3.86 ±0.81) X 10"^ 


4.781 


Cll 


2.17 X 10-3 




C17 


(3.93 ±0.78) X 10-3 


5.031 


C21 


2.47 X 10-3 




C22 


3.31 X 10-3 




C28 


1.28 X 10-3 




C32 


1.52 X 10-3 




C38 


(2.28 ±0.65) X 10-3 


3.503 


C40 


1.23 X 10-3 




C42^ 


2.47 X 10-3 




C44 


1.08 X 10-3 




C45 


3.19 X 10-3 




C46 


1.08 X 10-3 




C47 


1.39 X 10-3 




C50 


1.72 X 10-3 




C5is 


1.19 X 10-3 




C52^ 


1.09 X 10-3 




C53 


1.06 X 10-3 




C54'5 


2.44 X 10-3 




C55^ 


1.43 X 10-3 




C56 


1.64 X 10-3 




C59 


1.80 X 10-3 




C60 


1.41 X 10-3 




C62^ 


2.30 X 10-3 




C63 


2.82 X 10-3 




C64^ 


3.51 X 10-3 




C65 


1.50 X 10-3 




C66 


1.95 X 10-3 




C67 


3.27 X 10-3 




C68 


1.71 X 10-3 




C69 


2.34 X 10-3 




C77 


(1.96 ±0.61) X 10-3 


3.187 


C84^ 


1.26 X 10-3 




C92^ 


1.96 X 10-3 




C93^ 


1.74 X 10-3 




C96 


2.43 X 10-3 




C97 


1.66 X 10-3 




C98 


(3.08 ±0.79) X 10-3 


3.893 


C99 


3.28 X 10-3 




ClOl 


1.74 X 10-3 




C106^ 


1.46 X 10-3 




C107 


2.58 X 10-3 




C108 


2.29 X 10-3 











Table 3: Detection of X-ray source in NGC 2516: the label or ^ indicates variable type according to the 
classification reported in the text, count rate or upper limits are expressed in count s-^. 
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the cluster have been observed. 

The new results confirm the normal incidence of 
shorter period variables (likely d Scuti) and the anoma- 
lous incidence of longer period variables previously 
found. The results have been also compared to the 
survey in the X band obtained by ROSAT but no re- 
lation between variables and X-ray sources could be 
found. 

The statistic used is designed to be less sensitive to 
spurious effects such as air transparency variations. 
Observational effects can therefore be excluded and 
the variability ascribed to variations of the luminosity 
of the sources. 

While the hypothesis of binarity for explaining 
longer period variables can not be ruled out through 
our data, the position occupied by them suggests their 
possible belonging to the newly discovered class of 7 
Dor variables. 

If confirmed, the 7 Dor nature of the longer period 
variables in NGC 2516 would be an important result 
because it is unique among clusters, and the 8 objects 
disentangled in our survey certainly deserve further 
investigation. 
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This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 



